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ABSTRACT 

The Gaia satellite is a high-precision astrometry, photometry and spectroscopic ESA cornerstone mission, cur- 
rently scheduled for launch in 2012. Its primary science drivers are the composition, formation and evolution 
of the Galaxy. Gaia will achieve its unprecedented positional accuracy requirements with detailed calibration 
and correction for radiation damage. At L2, protons cause displacement damage in the silicon of CCDs. The 

resulting traps capture and emit electrons from passing charge packets in the CCD pixel, distorting the image 
PSF and biasing its centroid. Microscopic models of Gala's CCDs arc being developed to simulate this effect. 
The key to calculating the probability of an electron being captured by a trap is the 3D electron density within 
each CCD pixel. However, this has not been physically modelled for the Gaia CCD pixels. In Seabroke, Holland 
& Cropper (2008), the first paper of this series, we motivated the need for such specialised 3D device modelling 
and outlined how its future results will fit into Gala's overall radiation calibration strategy. In this paper, the 
second of the scries, we present our first results using Silvaco's physics-based, engineering software: the ATLAS 
device simulation framework. Inputting a doping profile, pixel geometry and materials into ATLAS and compar- 
ing the results to other simulations reveals that ATLAS has a free parameter, fixed oxide charge, that needs to 
be calibrated. ATLAS is successfully benchmarked against other simulations and measTirements of a test device, 
identifying how to use it to model Gaia pixels and highlighting the effect of different doping approximations. 
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1. INTRODUCTION 

The primary science drivers of Gaia, the ESA cornerstone mission scheduled for launch in 2012, are the com- 
position, formation and evolution of the Galaxy. Gaia will map the positions and kinematics of 10^ Galactic 
stars complete to F = 20 mag with unprecedented precision. However, Gala's required centroiding accuracy per 
observation, on the order of one-thousandth to one-tenth of a pixel dependent on magnitude, is already chal- 
lenging the calibration of today's state-of-the-art CCDs. Gala's CCDs were specifically designed for Gaia and 
are currently being manufactured by e2v technologies (UK). There are many additional calibration challenges 
facing Gaia^ (see Ref. 1 for a more detailed discussion of Gala's operational aspcc;ts and tests of Flight Model 
CCDs with respect to calibration needs). This paper focuses on the calibration challenge of needing to correct 
for radiation damage. 

Operating at L2, the satellite is expected to experience an end-of-mission (after five years) particle fiuence 

of X 10^ (mainly solar) protons cm~^ (10 MeV equivalence). In the silicon (Si) of the 106 Gaia CCDs, 
this will cause displacement damage in the form of electron traps. The CCDs are operated in Time Delay and 
Integrate (TBI) mode because the satellite will continuously spin (and precess). Depending on trap occupancy 
level determined by illumination history, traps may capture electrons from charge packets in pixels in the leading 
edge of image point spread functions (PSFs). Slow-emission (compared to the TDI period) traps can release 
their captured electrons after the image has passed. This charge loss from the image leads to lower image signal- 
to-noise. Periodic charge injection can keep the slow traps full but cannot fill the fast-emission traps^ (see Ref. 
2 for more details of the trap types). Fast traps can release their captured electrons into charge packets in pixels 
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in the PSF trailing edge. This will distort the PSF shape such that it no longer matches the fitting function 
used to centroid the image. The distorted PSF leads to a shifted centroid. This effect has been measured to be 
up to on the order of one-tenth of a pixel, which, depending on the brightness of the object, can be greater than 
the required centroiding accuracy. 

Gaia's approach is to have a radiation damage correction model that can produce a range of distorted fitting 
functions caused by a range of damage levels. The current approach does not fit a single distorted PSF to the data 
but instead must produce a damaged image at each iteration of the image parameter determination process. This 
is because the damaged profile depends on the sub-pixel location of the source and so is tied closely to the 'true' 
image parameters. The most accurate radiation damage model is assumed to be a microscopic implementation 
of the Shockley- Read-Hall theory^' ^ that models individual traps and charge packets. Monte Carlo techniques 
are required to model the stochastic probability of traps capturing and releasing electrons and to average over a 
range of trap positions within a pixel. The probability of an empty trap capturing an electron depends on the 
electron velocity (determined by known fundamental constants), trap cross-section (approximately known for 
some traps but often left as a free fitting parameter) and electron density in the vicinity of the trap (unknown 
but charge packets within pixels can be physically modelled by simultaneously solving for electrostatic potential 
and charge density using the Poisson and charge continuity equations respectively). Ongoing radiation test 
campaigns irradiate Gaia CCDs providing limited parameter space test data from which the effects of radiation 
damage on Gaia data have been derived and with which microscopic models can be verified. Subsequently, 
microscopic models can provide full parameter space synthetic data to test Gaia's data processing pipeline. 

However, microscopic models are currently too slow to be included in Gaia's data processing pipeline due 
to the extremely high data aquisition rate. Instead, much faster and consequently simpler macroscopic (e.g. 
column integrated) models are required to conduct the radiation damage correction to the data. These have 
developed from microscopic models and require microscopic models (as well as test data) to rigorously test 
them. One of the proposed macroscopic models^ is based on how the volume of a charge packet grows with its 
number of constituent electrons. Hence, both microscopic and macroscopic models require knowledge of electron 
distribution as a function of position within a pixel and as a function of the number of constituent electrons as 
inputs. Therefore the modelling of these inputs for Gaia's e2v CCD91-72 pixels using commercially available, 
physics-based, three dimensional (3D) device simulation software by Silvaco was proposed.^ 

This paper, the second in the series, presents the first results: benchmarking the Silvaco software. Because it 
is much easier to physically measure electrostatic potential (</>) than electron distributions, benchmarking of this 
type generally consists of simulating a device's 4> distribution and comparing it to that measured in the device. 
As (p determines the electron distribution, if the simulated (j) agrees with reality, this suggests the simulated 
electron distribution will also agree with reality. Silvaco is benchmarked against other simulation software and 
physical measurements. Section [2] introduces e2v's process parameters, which guide the benchmarking exercise. 
The results are presented in Sec. |3] with conclusions and future work in Sec. |4j 

2. E2V PROCESS PARAMETERS 

If the voltage difference between an electrode or gate (g) and the substrate (ss) is Vgss and the voltage on 
the underlying p-type substrate is Vss, then the cj) distribution through the buried channel (BC) part of the 
structure is as shown in Fig. [ij The maximum in the (buried) channel is designated 4>q]^- In the case of no 
signal charge being stored: (p^-^ ~ SVgss + ^cJiO' "^^^re (f>chO ^^"^ value of 0^,^^ for Vgss = volts (V). (p^-^ is 
termed the "channel parameter" and essentially relates the voltage in the underlying channel to that applied to 
the electrode. The actual value depends on the oxide capacitance per unit area, the doping distribution in the 
channel, the substrate doping concentrations and any fixed charge in the oxide. The factor S also depends on 
these quantities, being unity for lightly doped substrates and reducing slightly as the doping increases. 

If the voltage on the electrode is taken increasingly negative of the substrate, there comes a point when the 
voltage at the insulator-semiconductor (silicon dioxide-Si) interface becomes equal to that of the substrate, which 
allows electron holes to flow in from the substrate and accumulate at the interface. The resulting distribution is 
as shown in the left schematic of Fig. [l] If the electrode voltage is taken further negative, more holes accumulate 
and 4>f^Yi does not change and is said to be "pinned" . This maximum pinned cj) value is designated (jjchv ^'^^ 
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Figure 1. Left: Schematic (j) profile through an insulating oxide layer (gate dielectric), buried channel (BC) and substrate 
when a CCD pixel is in Inverted Mode Operation (IMG). The line between the n-type BC and p-type substrate is the p-n 
junction. Right: Silvaco TonyPlot comparing realistic and step (box) doping profiles: net dopant density as a function 
of depth into the device. The background is colour-coded according to the material structure of the ATLAS models as 
a function of depth: conductor (polysilicon electrodes, depth 0-0.5 ^m), Si02 (insulator, depth 0.5-0.63 nm) and silicon 
(semiconductor, depth 0.63-16 fim). Because the doping density is absolute, the minimum point on the realistic doping 
profile is the p-n junction, separating p-type ions to the right and n-type ions to the left. 

the gate voltage at the onset of pinning is designated Vgssp- The quantities ^^JjO' "^chp ^'^^ ^ssp can all be 
measured using test transistors included around the periphery of a device and values thereby provide useful data 
for process control (and benchmarking device simulations). 

3. BENCHMARKING 

e2v fabricated a test structure, a custom field effect transistor (FET), called TSlOO from a CCD66 (non-Gaia) 
wafer. Figure |2] shows that TSlOO consists of 16 transistors in parallel, sharing a common source but with 16 
different widths of gate G2 across cross-section A-A, leading to 16 individual drains. The experimental setup 
involved applying Vqi — —10 V and Vq2 = Vss = V. The channel created by the apphed gate voltages is 
essentially the same as a BC in a CCD, where electrons flow along ^^j^ from the common source to each drain. 
Measuring the onset of current flow gives '/'(.jj- 

A simulation of TSlOO should be able to reproduce the measurements of ^ function of the width (W) 

of gate G2: the presence of Gl increasingly reduces ^ the width of G2 across cross-section A-A reduces. 
The range of W covers the range of pixel feature sizes in the Gaia CCDs: in the Gaia along (AL) scan (charge 
transfer) direction, the electrode widths are 2 and 3 /im but the potential below the gates is always formed by 
two adjacent gates (i.e. combined width of 5 fim due to the simultaneous clocking scheme); in the Gaia across 
(AC) scan direction, the anti-blooming drain (ABD) is 2 /im wide, the supplementary BC doping (SBC) widths 
are 3 and 4 /im, the ABD shielding is 4 /im wide and the BC doping is 21 /im wide.^'^ 

The modelling input that electron distribution and (p are most sensitive to is doping. e2v modelled their 
doping implant process using SSuprem3, the one-dimensional (ID) process simulation module in ATHENA, 
which is Silvaco's process simulation framework that "enables process and integration engineers to develop and 
optimize semiconductor manufacturing processes" . This yielded a spread of simulated doping proflles: net dopant 
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Figure 2. Left: Microscopic image of TSlOO on a CCD66 wafer. The connection to the common source is labelled SRC (just 
above centre of the image) and common source itself runs the entire length of the transistor. SS labels the connection 
to the underlying substrate made on the front surface of the device. The overlying gate (labelled BG for bias gate and 
annotated Gl in the schematic) is visible along the entire length of the transistor. The increasing width of the underlying 
gate (labelled CHG for channel gate and annotated G2 in the schematic) from left to right is only apparent due to the 
changing topography of Gl. The widths of the narrowest gates are barely visible at the left end of the transistor but the 
individual drain connections below the transistor indicate the gate positions. Each individual drain connection is labelled 
with the width of the gate. The drain connection to the narrowest gate is labelled 02 (2 fim) just to the left of BG above 
the transistor and connected to the drain below the transistor. The drain connection to the next narrowest gate is labelled 
02.5 just to the left of 02 above the transistor and connected to the drain below the transistor to the right of the 02 
drain. The labelling of the drain connections to gates 03 and 03.5 is off the left edge of the image above and below the 
transistor respectively. Drain connections to gates 04, 04.5, 05, 06, 07, 08, 09 and 10 are labelled along the bottom of 
the image. The labelling of the drain connections to gates 12, 14, 16 and 20 is off the right edge of the image. Right: 
Schematic of TSlOO, showing the narrowest 3 gates: 2, 2.5 and 3 pm. W is the width of gate G2 across cross-section A-A. 
A vertical view through cross-section A-A is at the top of the schematic where gate Gl overlays gate G2. 



density (p, donor ions cm~'^) as a function of depth (z, cm) into the device's Si layer. The resulting spread was 
calibrated against Secondary Ion Mass Spectrometry (SIMS) measurements and empirically fit by a function of 
the following form: 

p = / X 10"^'+"-+- + pss where Pss^^-j,, (1) 

where / is the dimensionless BC scale factor, a, b and c are constants proprietary to e2v, pss is the substrate 
doping density (acceptor ions cm^^), q is electron charge (C), pp is hole mobility (480 cm^V~^s~^) and R is the 
substrate resistivity (0cm). 

R = 20 ilcm for TSlOO, which gives pss = —6.51 x 10^^ acceptor ions cm"'^ (by convention donor and acceptor 
ions have positive and negative densities respectively because they are positively and negatively charged and thus 
can cancel out the effect of each other). Note the Gaia CCDs use much higher resistivity Si: R = 100 ficm 
for the Astrometric Field (AF) and Blue Photometer (BP) CCDs and R = 1500 ilcm for the Red Photometer 
(RP) and Radial Velocity Spectrograph (RVS) CCDs (A. Walker (e2v), private communication). Nevertheless, 
like the Gaia AF and BP CCDs, TSlOO also has the e2v standard depth of 16 /im, whereas the Gaia RP and 
RVS CCDs are e2v deep depletion devices with thicknesses of 40 /im. Although the simulations presented in this 
paper are for a lower resistivity device (TSlOO) than the Gaia CCDs, the same methods presented here will be 
used to simulate the Gaia CCDs and so the benchmarking of these methods also applies to simulations of the 
Gaia CCD pixels. 
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/ is a free parameter that cannot be easily physically measured. Instead, it is varied in Equation [T] to 
produce a doping profile, which is substituted into Poisson's equation to solve for (f). The only part of the 
simulated internal (p distribution that can be compared to measurement is <^(jjj. However, the difference between 
the measured and predicted 4>^^^ depends on flat-band voltage (^^), which cannot be independently measured. 
Therefore, models need to predict a parameter that is independent of Vq^ and can be measured: ^^jip- Using 
e2v's in-house simulation software, written by M. Robbins,^ / is varied to satisfy 4>q]^, ~ "^^lich is the 

standard value measured by e2v. This is assumed to apply to TSlOO because ^^hp ^'^^ strongly dependent on 
R. The resulting value of / is used in Equation [l] to derive a doping profile applicable to TSlOO, which is read 
into ATLAS, Silvaco's device simulation framework: "ATLAS is a physically based two and three dimensional 
device simulator. It predicts the electrical behaviour of specified semiconductor structures and provides insight 
into the internal physical mechanisms associated with device operation" . The primary user-inputs to ATLAS 
are device geometry, materials and the doping profile (s). 

Robbins' model implements the simplest possible numerical method to solve Poisson's equation in ID so it 
does not attempt to physically model V^. To derive requires modelling the test transistor that e2v used 
to measure Vgssp. Robbins models this in ID by applying voltage to a layer of Si02 on top of a slab of doped 
Si. Figure [3] shows that the Robbins model pins at Vgssp = —6.97 V, whereas e2v measured Vgssp ~ —8.5 V. 
This value was not measured from TSlOO itself but is an average value of the batch, which is consistent with 
the batch-to-batch average of —8.1 V for i? = 20 ficm. There is a spread of channel parameters wafer-to-wafer 
and batch-to-batch so the actual value that TSlOO actually pins at is not known. Nevertheless, this is the best 
information available with which to calibrate the Robbins model in order to approximately model TSlOO. If 
"Robbins realistic doping — 1.53 V" is plotted in Fig. |3]to represent a model that does include it shows 
that the Robbins model gives the same as the Robbins realistic doping line when T^ss — is applied to the 
gate, where = —1.53 V (by convention Vj^ is a negative number). ATLAS physically models Vj^ via the 
free parameter fixed oxide charge (Q), which is defined at all insulator-semiconductor interfaces per unit area (C 
cm^-^). 

The same Robbins test transistor model doping profile and oxide thickness was used in the corresponding 
ATLAS model, which has to be two-dimensional (2D). Therefore the main dimension is depth with an arbitrary 
distance in an orthogonal device dimension. The structure does not change in this direction and so the ATLAS 
simulation is planar. Figure [31 shows the models exhibit similar Non- Inverted Mode Operation (NIMO, T^ss > 
Vgssp) behaviour: the "ATLAS realistic doping: Q — 0" line has the same NIMO gradient as the Robbins model, 
only offset by «0.2 V. It shows that Q needs to be calibrated to the "Robbins realistic doping — 1.53 V" line in 
order for ATLAS to model = -1.53 V. 

The bottom plot of Fig. [3] highlights a difference in behaviour between the Robbins model and ATLAS in 
Inverted Mode Operation (IMO, Vgss < Vgssp). The Robbins model produces straight IMO lines because it is an 
electrostatic model, simulating </> differences due to space charge regions. ATLAS does not produce straight IMO 
lines because it is a full semiconductor model, simulating electron and hole carrier concentrations. The second 
order effect of more and more holes flowing to the Si-Si02 interface causes the ATLAS IMO lines to curve. This 
means it is more difficult to deflne ^(^jjp and Vgssp from ATLAS simulations. Because of this and the fact that 
the TSlOO G2 gates (and Gaia CCDs) are run in NIMO, Q is varied to align the NIMO parts of the ATLAS line 
and the "Robbins realistic doping - 1.53 V" line. Q = 2.26 x 10" C cm brings them into excellent agreement. 
However, the models cannot be compared in an absolute sense to the TSlOO measurements because TSlOO's 
T^ssp is not known, which means that TSlOO's cannot be derived for the Robbins model nor TSlOO's Q for 
ATLAS. Nevertheless, exactly calibrating Q to the Robbins model has the advantage of allowing a direct 2D 
comparison of ATLAS and the Robbins model, while only being able to compare to TSlOO measurements in a 
relative sense. 

The average doping profile in Equation [T] only applies to the BC part of the pixel and not profiles through 
SBCs and ABDs, which are features in the Gaia image pixels. To avoid doping and </> discontinuities in the 
AC pixel dimension, ideally ATLAS would use a 2D doping file derived from 2D ATHENA simulations as input 
(doping does not change in the AL direction and thus 3D doping would be superfluous). However, e2v have 
not modelled the implant processes for the BC, SBC and ABD simultaneously in SSuprem4, ATHENA's 2D 
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Figure 3. Plots of Robbins ID and ATLAS planar test transistor models simulating (p^^i against Vgss, showing ^^j^q, 
and Vgssp (top) and zoomed in on the ^^hp ^gssp region (bottom). 
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process simulation software. In the absence of such a 2D doping file, the easiest and quickest way to simulate 
complex doping profiles is using step profile doping. Realistic doping profiles can be approximated by uniformly 
doped boxes: e.g. the SBC part of the pixel would have one box for the SBC itself and one for the substrate 
(the constant density of which is given by Equation [ij , which is less than the SBC box doping forming a step 
profile. However, a BC reahstic doping profile adjacent to a SBC step doping profile would cause doping and </> 
discontinuities. Thus, the whole Gaia image pixel needs to be modelled with step doping. Both realistic and step 
doping ATLAS models are benchmarked in this paper to highlight the effect of different doping approximations. 

The TSlOO BC box parameters are the BC doping (n-type, implant donor dopant density Nj^, cm^'') and 
the BC implant depth {d, cm), over which N-q is uniformly distributed. ATHENA modelling by K. Ball (e2v) 
finds d = 0.54 yum. Like / in the realistic doping case, the value of N-q needs to satisfy '/'(.jjp = 4.6 V. Nj) is 
derived using the step doping model, which, assuming full depletion, allows the (f> distribution within the gate 
insulator and Si region to be analytically obtained.^ Using e2v's in-house simulation software implementation of 
the step doping model (written by M. Robbins), N-q — 2.76 x 10^^ cm~^. The right plot of Fig. ijcompares the 
resulting step doping profile with the realistic one. The BC box in the step doping model shoulohave the same 
number of donor dopant ions as the realistic doping profile. Integrating the BC part of Equation [l] (by excluding 
the pss term) with respect to z calculates the concentration of dopant ions in the BC per unit area of the device 
in the plane perpendicular to z (pgf^, cm~^). Nj) = pQ(^/d = 2.65 x 10"'^^ cm~'^, which is consistent with the 
uniformly doped BC box value. 

Figure |3] shows the step doping models also exhibit similar NIMO behaviour: the "ATLAS step doping: 
Q = 0" line has the same NIMO gradient as the Robbins step doping model. The different methods the two 
models employ, the Robbins step doping model is a ID analytical solution of Poisson's equation while the ATLAS 
planar step doping model is a 2D numerical solution of Poisson's equation, should produce the same results. The 
offset of «0.3 V shows fairly good agreement. With no Vq^ or Q, Vgssp(step) < Vgssp (realistic). Therefore the 
required to shift the Robbins step doping line to pin at —8.5 V is less than the realistic doping value of 
V^. The resulting "Robbins step doping — 0.77 V" line has a slightly steeper gradient than the realistic doping 
line. Because the different doping models are diverging around (IJqIjQ, the ATLAS step doping model has its Q 
calibrated against the "Robbins realistic doping — 1.53 V" to make the step doping model as similar as possible 
to the realistic doping model around which gives Q = 0.54 x 10^^ C cm^^. However, Fig. |3]shows that the 

excellent agreement between the ATLAS step and realistic doping models around (/'(^hO expense of good 

agreement around ^cj^p- The bottom plot of Fig. |3j shows that the IMO (f)^-^ of the ATLAS step doping model 
agrees with the Robbins models at ^4.6 V, whereas the realistic doping models are offset by —0.1 V at ^4.5 V. 

"The reference (j) can be defined in various ways. For ATLAS, this is always the intrinsic Fermi (p {(j)^)" (ATLAS 
Users Manual). In an intrinsic semiconductor, 0j — (j)^, where (j)^ is the Fermi (j), which is the electrochemical (j) 
equivalent to voltage. However, in extrinsic (doped) semiconductors, referencing (p^ to 0j means that = V 
in the following equation, which only applies in a p-type substrate: 



where k is Boltzmann's constant (J K^^), T is temperature (it is assumed that TSlOO measurements were 
conducted at room temperature: 293 K), Np^ is acceptor doping density {Np^ — \pss\) ^^nd Tij is the intrinsic 
carrier density («7 xlO® cm^^). Substituting the values for TSlOO into Equation [2] gives <f)^ = —0.28 V, which is 
the same value seen in the substrate of the 0£ profiles in Fig. |4j With (/)ss = -0.28 V, for Vgss to equal zero, Vg 
must also equal —0.28 V. The Robbins model profiles have their 4> referenced to (p^s — 0. Therefore, to directly 
compare both model outputs, the ATLAS simulations were also referenced to 0ss = by adding 0.28 V to every 
mesh point through the gate, oxide and Si of the ATLAS profiles and plotted in Fig. |4] (this was also done for 
Vgss and i/i^,-^ for every ATLAS simulation in Fig. |3]). 

The bottom plot of Fig. l4| shows that the distribution around (j)^Yi differs between step and realistic doping. 
Both the step doping models produce a shallower and narrower peak in terms of depth into the Si. The Robbins 
step doping model produces a higher than both the realistic doping models. Because the step doping (j) 
distribution around (j)^Yi is constant in shape and depth and is scaled in by Q, it is possible to remove one of 
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Figure 4. Plots of Robbins ID and ATLAS planar test transistor models simulating NIMO against depth into device, 
comparing the resulting from realistic and step doping profiles (top) and zoomed in on the f/igjjQ region (bottom). The 
top plot shows the full depth of TSlOO: electrodes (depth 0-0.5 /im), Si02 (depth 0.5-0.63 ^m) and Si (depth 0.63-16 
/im). The Si-Si02 interface is 0.63 /im deep. 
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the differences between the step and reaKstic doping distributions by calibrating Q to mafce the ATLAS step 
doping model 0^,^^ agree with the (t)f,^^ of the realistic doping models. The other differences in the 4> distributions 
around 0^,^ are intrinsic to the differences between the realistic and step doping profiles. 
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Figure 5. Silvaco TonyPlots of the 2D ATLAS model showing the two ends of TSlOO: middle gate G2 with a width of 
2 ^m {left) and 20 /im [right) with (f)^ contours overlaying the electrodes, Si02 layer and the Si, due to Vq-^ — —10.28 V 
either side of Vq2 ~ —0.28 V. The gaps between the electrodes are 0.25 fim and filled with oxide. 

Having derived the realistic and step doping profiles and from f D Robbins simulations and Q from planar 
ATLAS simulations, these values can now be applied in 2D models of TSf 00. TSf 00 only needs to be modelled 
in 2D because the length of gate Gf in Fig. [2] in the charge flow direction is long enough for 3D effects to be 
negligible. Both the 2D Robbins model^ and ATLAS model simulate TSf 00 by modelling the cross-section A-A 
in Fig. |2]but only including one G2 gate width at a time. The cross-section A-A shows that Gf overlaps G2. 
This does not need to be modelled because only conductors in contact with the insulator over the Si determine 
the 4> distributions in the Si. Therefore, Fig. |5] shows how the gates were simulated and how the (j> under the Gl 
gates pinches off the </> under G2 when G2 is narrow (left plot) compared to when it is wide (right plot), (p^^i 
is measured from a (f> profile through the centre of each G2 gate width and is plotted in Fig. [6] fn ATLAS, the 
maximum (/> under G2 is actually 0^,^^ — 0ss, where 0ss — —0.28 V. 

The models simulate (/)(internal), which cannot be directly measured. 0(transistor measured) is less than 
(/)(internal) because it does not take into account the "diode drop" . This is the 4> required to surmount the p-n 
junction barrier (diode) and commence conduction, called built-in (jj {(f'h])- Therefore, (^(measured) does not 
measure the diode drop because it has "used up" (j)^^ to allow the measurement to be made and so (j)^^ needs 
to be added back on to the (/)(measured) to make it up to the value of ^(internal), to compare directly with 
the models. The exact value of cj)-^^ depends on the doping profile of the transistor's source and drain but the 
standard value of 0.60 V is assumed. Thus all the e2v measurements in Fig. [6] have had 0.60 V added to them. 

Figure [61 shows all the models and measurements have (j)^-^ values that plateau at large G2 gate widths. The 
top plot oiFig.[6]shows that at W = 20 /im, the ATLAS realistic and step doping 2D models give 0^,^ = f f .50 V. 
This is the same as the (PqI^q — f f .50 V produced by the planar ATLAS (and fD Robbins) models in Figs. |3] 
and |4j However, one may have expected a difference in between the ATLAS planar and 2D simulations 
because the planar simulation does not include any fringing fields (because there is only one gate), whereas the 
2D simulation does include fringing fields, which can decrease (j^di- This decrease is seen comparing the Robbins 
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Figure 6. Plots of (i)^ as a function of G2 gate width, comparing measurements of TSlOO by e2v to 2D simulations of 
TSlOO by the Robbins and ATLAS models, absolutely (top) and relatively [bottom). 
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ID and 2D models: the top plot of Fig. |6]gives (p^^i = 11-26 V at W = 20 iim, which is 0.24 V less than the (^^hO 
= 11.50 V produced by the Robbins ID model. The difference between the ATLAS and Robbins 2D models 
could be related to how the models simulate (p in the oxide gaps between the electrodes: ATLAS solves Poisson's 
equation at every mesh point, whereas the Robbins model linearly interpolates between the electrodes. The 
ATLAS realistic doping model produces rvO.l V less than the Robbins models in ID at Vgss — —10 V in Fig.jSj 
suggesting in 2D one may have expected the Robbins model to produce higher cj) values than ATLAS, rather 
than vice versa as seen in Fig. [6j 

The ATLAS step doping model, which was deliberately calibrated to the ATLAS planar realistic doping model, 
shows excellent agreement with the ATLAS realistic doping model in 2D at large gate widths but diverges as 
W decreases (and fringing fields become closer in size to the gate width) until reaching its maximum offset of 
^0.2 V at W = 2 fim. This suggests that ATLAS models using step doping profiles do not simulate fringing 
fields as well as models using realistic doping profiles. Gaia charge packets will sit in (j> defined by two adjacent 
gates with combined width of 5 /im wide so modelling the Gaia pixel with step doping could introduce a (j) offset 
of 0.1 V but this is much less than the expected variation of up to ±15% from Gaia CCD production spreads.® 

The absolute differences between the models are not as important as the relative differences because electrons 
(and holes) only "see" relative changes in cj) rather than absolute values. The bottom plot of Fig. [6] shows there is 
excellent agreement between all the models at large gate widths when each model is normalised to its (/)^^ when 
W = 20 /im. The Robbins and ATLAS models diverge at smaller gate widths to a maximum offset of wO.3 V 
at W = 2 /im. Again, this divergence is indicative of the differences of how the models simulate the oxide gaps 
and the 2D effect of fringing fields. 

The gate widths plotted in Fig. |6] are the design values. However, manufacturing tolerances mean that all 
dimensions of the electrodes can be larger or, more likely, smaller than the nominal value by up to 1 /im due to 
under- and over-etching respectively. All gates will differ from their nominal dimensions by the same amount. 
Plotting the e2v measurements minus 1 /im in the bottom plot of Fig. [6] brings the measurements into much 
better agreement with both the Robbins and ATLAS realistic doping models than previously and provides a 
plausible explanation for the relative 4> offset between the measurements and the models. 

4. CONCLUSIONS AND FUTURE WORK 

The aim of this series of papers is to derive the electron distribution within a charge packet, as a function of 3D 
position within the different Gaia CCD pixels and as a function of the number of electrons in the charge packet, 
to feed into Gaia radiation damage correction models. This paper has presented the first results of using the 
Silvaco ATLAS device simulation software, benchmarking it against other simulation software (Robbins' models) 
and physical measurements of a test structure. 

Inputting the correct doping profile, pixel geometry and materials into ATLAS and comparing to the corre- 
sponding ID Robbins model reveals that ATLAS has a free parameter, fixed oxide charge {Q) that needs to be 
calibrated. Because the measurement required to calibrate ID models to e2v's test structure TSlOO (Vgssp) was 
not made, Q was calibrated against the Robbins ID model. This was done in the Non-Inverted Mode Operation 
(NIMO) part of the simulated <p(.]j-Vgss curve because of the different behaviour of the models in the Inverted 
Mode Operation (IMO) part and the fact that the TSlOO G2 gate and Gaia CCDs operate in NIMO. 

Comparison of ATLAS and Robbins model potential {(p) profiles highlights the different reference potentials 
{(p) of the models. ATLAS is referenced to the intrinsic Fermi cj), while the Robbins model is referenced to 
(/-ss = 0. For TSlOO with Vgss(Gl) = -10 V, %s(G2) = V and V^ss = V, this meant the ATLAS model 
apphed ^g(Gl) = -10.28 V and Vg{G2) = -0.28 V because (jjss = -0.28 V even with Fss = V applied. Hence 
to simulate the Gaia 4> distribution under Vg(high) = 10 V and Vg(low) = V (Vss = V) will require applying 
Vg + (pss to each gate. This will mean that electrons in the Gaia 3D pixel ATLAS simulations will sit in slightly 
different absolute values compared to reality. The electron distribution will not differ from reality as long as 
the relative (j) values between adjacent electrodes are consistent. 

Because ATLAS was calibrated to the Robbins model in ID, this permitted a direct absolute comparison of 
the models in 2D. This reveals a small (p^.^^ difference between the models in 2D, ranging from 0.24-0.5 V going 
from the widest to the narrowest gate. As the Vj^ and Q used in the 2D models are not specific to TSlOO, the 
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Robbins and ATLAS models could only be compared to measurements in a relative sense. While the models 
differ slightly in their predicted (^^j^ under narrow gate widths, the general trend of both models follows the 
physical measurements of TSlOO well (after correction for a suspected systematic offset in TSlOO itself). 

It was found that step doping produces different (f) distributions around (j)^^^ to realistic doping. As the 
realistic doping is more accurate than step doping, the produced by step doping will not be as physically 
realistic as if realistic doping had been used. Consequently, the electron distributions produced by step doping 
may not be as physically realistic as if realistic doping had been used. However, it is not yet clear to what extent 
changes in (p affect electron distributions. This is will be investigated in the course of this work. In the absence of 
realistic doping profiles for Gaia image pixels, these pixels will be modelled with step doping profiles. Although, 
this will provide only an approximation of physically-derived electron distributions, it will be an improvement 
on arbitrarily-chosen electron distribution functions used in previous Gaia radiation damage models.^ 

Therefore, ATLAS has been successfully benchmarked against a device with lower resistivity (TSlOO) than 
the Gaia CCDs. Because the simulation methods arc the same, regardless of resistivity, the benchmarking results 
are equally applicable to Gaia simulations, identifying how to use ATLAS to model Gaia pixels and the effect of 
using step rather than realistic doping. This is the first attempt by anyone to model the Gaia pixel to derive its 
electron distributions. It may also be the first time this has been done for any e2v pixel. If this is not the first 
attempt by anyone to model an e2v pixel, it is the most complex e2v pixel to be modelled, as the Gaia image 
pixel is the most complex pixel architecture ever made by e2v, because it includes so many different features: 
buried channel, supplementary buried channel and anti-blooming drain and shield. A realistic microscopic Gaia 
pixel model will not only provide inputs to Gaia radiation damage correction models but be able to simulate 
and physically explain the results found by detailed measurements of close- reject Gaia Flight Model CCDs.^ 
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